Method and apparatus for distributing a quantum key

ABSTRACT

A method for distributing a quantum key is provided, including sending a first photon to a first receiver; sending a second photon to a second receiver, the first and second photons being a pair of time-energy entangled photons; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins, wherein a combination of a time bin and a frequency bin corresponds to a character.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority to U.S. provisional application 61/543,691 filed Oct. 5, 2011, which is incorporated by reference in its entirety herein.

FIELD

The disclosed subject matter generally relates to integrated optical systems. More particularly, the disclosed subject matter relates to a new design for a photonic crystal useful for a wide variety of applications.

BACKGROUND

Quantum key distribution (QKD) enables two parties to establish a secure key at a distance, even in the presence of one or multiple eavesdroppers. The key can then be used to secretly transmit information using unconditionally secure one-time pad encryption. Unconditional security is achieved by the laws of physics rather than by assumptions about the computational abilities of the eavesdropper. Various QKD protocols have been proposed and implemented. A first protocol due employed polarization states of photons passed between the two parties. There has been growing interest in schemes employing photons in Hilbert spaces of high dimension, resulting in a potentially very large alphabet size. Different degrees of freedom have been considered, including polarization, time, and spatial modes.

What is needed is a protocol that allows two parties to generate their secure key at the maximum rate allowed using the time-energy basis. Moreover, a protocol is needed that is compatible with fiber-based dense WDM (DWDM) systems commonly used in classical fiber communications. It is also desirable to have an extremely compact, stable, and scalable platform for the protocol's implementation.

SUMMARY

A method for distributing a quantum key is provided, including sending a first photon to a first receiver; sending a second photon to a second receiver, the first and second photons being a pair of time-energy entangled photons; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins, wherein a combination of a time bin and a frequency bin corresponds to a character.

In some embodiments, the first photon and the second photon are sent via optical fiber. In some embodiments, the first photon and the second photon are sent via a photonic integrated chip. In some embodiments, the first photon and the second photon are sent through free space.

A method for distributing a quantum key is provided including generating a plurality of pairs of time-energy entangled photons; for each pair of time-energy entangled photons, sending one photon to the first receiver and one photon to the second receiver; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins wherein each pair of one of the plurality of time bins and one of the plurality of frequency bins corresponds to one of a plurality of characters.

A method of receiving a quantum key is provided including receiving a first photon, the first photon being one of a pair of time-energy entangled photons; detecting a time of arrival of the first photon; detecting a frequency of the first photon; and determining a character based on the detected time and frequency of the first photon.

In some embodiments, the determining a character step further comprises: assigning the detected time to a time bin; assigning the detected frequency to a frequency bin; and determining the character from a coding scheme based on the time bin and the frequency bin.

In some embodiments, the method further includes detecting whether an eavesdropper has detected the first photon by determining the frequency distribution of the first photon.

A method of receiving a quantum key is provided including receiving a plurality of photons, each of the plurality of photons being one of a pair of time-energy entangled photons; detecting the time of arrival of each of the plurality of photons; detecting the frequency of each of the plurality of photons; and determining a character based on the detected time and frequency of each of the plurality of photons.

An apparatus for receiving a quantum key is provided including a channel configured to receive photons; a plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; and a plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to send an indication of photon detection.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1( a) illustrates a system using the time-coding scheme in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( b) illustrates time bins in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( c) illustrates a further system using the time-coding scheme in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( d) illustrates time bins in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 2( a) illustrates Franson interferometers in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 2( b) illustrates the coincidence counting probability P_(C) displayed versus δt in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 3( a) illustrates mutual information as a function of detector jitter, σ_(det), normalized to σ_(bin) in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 3( b) illustrates mutual information as a function of frequency channels, n_(f) in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 4 illustrates a further Franson interferometer in accordance with exemplary embodiments of the disclosed subject matter.

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS

It is understood that the subject matter described herein is not limited to particular embodiments described, as such may, of course, vary. It is also understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present subject matter is limited only by the appended claims. Where a range of values is provided, it is understood that each intervening value between the upper and lower limit of that range, and any other stated or intervening value in that stated range, is encompassed within the disclosed subject matter.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosed subject matter belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present disclosed subject matter, this disclosure may specifically mention certain exemplary methods and materials.

All publications mentioned in this disclosure are, unless otherwise specified, incorporated by reference herein for all purposes, including, without limitation, to disclose and describe the methods and/or materials in connection with which the publications are cited.

The publications discussed herein are provided solely for their disclosure prior to the filing date of the present application. Nothing herein is to be construed as an admission that the present disclosed subject matter is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided may be different from the actual publication dates, which may need to be independently confirmed.

As used herein and in the appended claims, the singular forms “a,” “an” and “the” include plural referents unless the context clearly dictates otherwise.

Nothing contained in the Abstract or the Summary should be understood as limiting the scope of the disclosure. The Abstract and the Summary are provided for bibliographic and convenience purposes and due to their formats and purposes should not be considered comprehensive.

As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosed subject matter. Any recited method can be carried out in the order of events recited, or in any other order that is logically possible.

Reference to a singular item includes the possibility that there are plural of the same item present. When two or more items (for example, elements or processes) are referenced by an alternative “or,” this indicates that either could be present separately or any combination of them could be present together, except where the presence of one necessarily excludes the other or others.

As summarized above and as described in further detail below, in accordance with the various embodiments of the present invention, a large-alphabet protocol is described herein referred to as a ‘Wavelength Division Multiplexed Quantum Key Distribution’ (WDM-QKD) that employs time-energy entangled photon pairs. Given a certain photon flux n and channel bandwidth ΔΩ, WDM-QKD allows two parties, e.g., “Alice and Bob,” to generate their secure key at the maximum rate allowed using the time-energy basis. WDM-QKD forms the quantum-cryptography analog of present-day WDM systems and is compatible with fiber-based dense WDM (DWDM) systems commonly used in classical fiber communications. An implementation of WDM-QKD is described herein as both a modern fiber optic network and in a photonic integrated chip (PIC). The latter provides an extremely compact, stable, and scalable platform for the protocol's implementation.

WDM-QKD enables Alice and Bob to generate their shared key at a maximum rate of n log₂(ΔΩ/n) at log₂(ΔΩ/n) bits per photon (bpp), using present-day single photon counters, WDM equipment, and simple components such as beam splitters or directional couplers. The information per photon per bandwidth is log₂ (ΔΩ/n)/ΔΩ.

QKD protocols derive security from the fact that a measurement changes an unknown state if measured in a conjugate basis. If Alice and Bob employ two conjugate bases, they can detect a measurement by a third party eavesdropper, “Eve,” including quantum nondemolition measurements. In the protocol described herein, the conjugate bases are time and energy. Alice and Bob make time and energy measurements using an extended Franson interferometer. Assuming 400 wavelength channels with detectors having 40 ps jitter, it is expected that for a bandwidth ΔΩ=10 nm around 1550 nm, Alice and Bob can generate a key at a maximum of 10 Tera bits per second (Tbps) at 1 bit per photon (bpp); or at 10 bpp and 100 Gbps; or at 200 Mbps at 20 bpp.

Time-energy entangled photon pairs with correlation time σ_(cor) generated by spontaneous parametric down conversion (SPDC) were considered herein, assuming a pump field at frequency ω_(p) with coherence time σ_(coh). The photon pair wave function can be written as |Ψ

=∫_(−∞) ^(∞)∫_(−∞) ^(∞)ψ(t _(A) ,t _(B))|t _(A) ,t _(B) ,o _(A) ,o _(B)

dt _(A) dt _(B)  (1) |t _(A) ,t _(B) ,o _(A) ,o _(B)

={circumflex over (α)}_(o) _(A) ^(†)(t _(A)(t _(A)){circumflex over (α)}_(o) _(B) ^(†)(t _(B))|0

ψ(t _(A) ,t _(B))∝e ^(−(t) ^(A) ^(−t) ^(B) ⁾ ² ^(/4σ) ^(car) ² e ^(−(t) ^(A) ^(+t) ^(B) ⁾ ² ^(/16σ) ^(car) ² e ^(iw) ^(p) ^(/2(t) ^(A) ^(+t) ^(B) ⁾

The creation operator {circumflex over (α)}_(o) _(i) ^(†)(t_(j)) denotes creation in spatial mode o_(i) at time t_(j). If the correlation time of the photons is sufficiently smaller than the time bin duration, Alice and Bob can establish a secret key using the correlated photon arrival time. A schematic of this protocol is shown in FIG. 1. The continuous biphoton wave function can be discretized as a sum over time bins |σ_(bin) ^(i)

of duration σ_(bin) as

$\begin{matrix} {{\left. \overset{\_}{\Psi} \right\rangle = {\sum\limits_{i = {- \infty}}^{\infty}{\sum\limits_{j = {- \infty}}^{\infty}{G^{i,j}\left. {\sigma_{bin}^{i},\sigma_{bin}^{j}} \right\rangle}}}}{G^{i,j,} = {\int_{i\;\sigma_{bin}}^{{({i + 1})}\sigma_{bin}}{\int_{j\;\sigma_{{bi}\; n}}^{{({j + 1})}\sigma_{{bi}\; n}}{{\psi\left( {t_{A},t_{B}} \right)}{\mathbb{d}t_{A}}{{\mathbb{d}t_{B}}.}}}}}} & (2) \end{matrix}$

The probability that Alice and Bob project into time bins σ_(bin) ^(i) and σ_(bin) ^(j) is therefore p^(i,j)=

σ_(bin) ^(i), σ_(bin) ^(j)| ψ

|²=|G^(i,j)|². Detector jitter of magnitude σ_(det) influences the fidelity of this projection, so that roughly σ_(bin)>σ_(det) for reliable communication. Jitter therefore reduces the maximum number of characters to σ_(coh)/σ_(det) from the maximum allowed, given by the Schmidt number, K≈σ_(coh)/σ_(cor). The fastest single photon detectors provide σ_(det)≈40 ps, whereas σ_(cor) can be on the order of 10 fs−1 ps.

FIG. 1( a) illustrates a system 100 using the time-coding scheme. A strong laser (not shown) pumps a nonlinear crystal. Photons pairs are generated by SPDC 102 and sent across channels of equal length to detectors 106 of Alice's computer 108 and Bob's computer 110 who measure their arrival times. FIG. 1( b) illustrates time bins agreed upon by Alice and Bob over a public channel. If a photon pair is detected in a given time bin, then that character is shared between Alice and Bob. FIG. 1( c) illustrates a system 200 in which photons pairs are also generated by SPDC 202. Alice's computer 208 and Bob's computer 210 include a grating, dispersive element or multi-channel filter 212 before their detectors 206 to obtain frequency information. FIG. 1( d) illustrates time bins.

For technological reasons, detector jitter is unlikely to approach the sub-ps regime in the near future. To overcome the technological mismatch between photon correlation time and detector jitter, the protocol described herein utilizes the circumstance that photons are not only entangled in time, but also in frequency. The biphoton wave function in the frequency domain, |Ψ_(F)

=FT₂|Ψ

, where FT₂ denotes the two-dimensional Fourier transform. Therefore, |Ψ_(F)

=∫∫

(ω_(A),ω_(B))|ω_(A),ω_(B) ,o _(A) ,o _(B)

dω _(A) dω _(B), where ψ(ω_(A),ω_(B))αexp[−σ_(cor) ²/4(ω_(A)−ω_(B))²]exp[−σ_(coh) ²/4(ω_(A)+ω_(B)−ω_(p))²].

States of this form show non-local frequency correlations. Thus, if Alice measures a frequency ω_(A) on one photon, then Bob must measure a frequency ω_(B)≈ω_(p)−ω_(A) on the other, where ω_(B)=ω_(p)−ω_(A) for σ_(coh)→∞. If both Alice and Bob place multi-channel filtering elements before their detectors (as shown in FIG. 1( c)), then they can collect frequency information in addition to timing information. This filtering can be modeled as Gaussian projections onto discrete output spatial modes |ζ^(i)

with frequency bandwidth δv and center frequency v_(i). Such filtering increases the correlation time of the photon pair, but this can be kept lower than the detector jitter to minimize the loss of timing information. Alternatively, one can increase the time bin duration σ_(bin) to allow for smaller δv. Optimization of these parameters given limitations of the source and transmission channel are described below.

The QKD protocol is secure as discussed herein. Measurements of the frequency and creation time of the biphoton packet disturb its wave function in these bases; temporal positive operator valued measurements (POVM) reduce the coherence time while frequency POVMs increase the correlation time as discussed in greater detail below.

Alice and Bob can check the security of their channel by testing for changes in the correlation and coherence times. They use an extended Franson interferometer (eFI) 300/301, depicted in FIG. 2( a), in which photons pairs are also generated by SPDC 302. eFI 300/301 include, e.g., beam splitters BS 326, partial beam splitters PBS 328, and detectors 306. Measurement in Alice's and Bob's arm of the eFI are described by the annihilation operators for the long and short paths, with respective times t_(L) ^(i) and t_(S) ^(i), and iεA, B denoting Alice and Bob: âA(tA)=1/√{square root over (2)}[âA(t_(L) ^(A))+âA(t_(S) ^(A))] and âB(tB)=1/√{square root over (2)}[âB(t_(L) ^(B))+âA(t_(S) ^(B))].

FIG. 2( a) illustrates for Bob's eFI 301, the short arm is switched between position A and B to achieve lengths δt₁ and δt₂. This allows determination of both and σ_(cor) so that weak frequency and time measurements on the photon pair can be detected. FIG. 2( b) illustrates the coincidence counting probability P_(C) displayed versus δt. For example, at position A, δt=0 and P_(c)=1.0.

These times are redefined using standard notation, where Δt=(t_(L) ^(A)−t_(S) ^(A)) and δt=(t_(L) ^(A)−t_(S) ^(A))−(t_(L) ^(B)−t_(S) ^(B)). A selection is made to use Δt, Δt−δt>>σ_(cor) in order to avoid single photon interference between long and short paths of a single arm of the eFI. Selection of two different values, δt_(1,2), is allowed by placing a switch in the long path of one arm of the eFI. The parameter δt_(1,2) is then scanned to provide phase-dependent coincidence counting measurements. Alice and Bob measure the coincidence counting probability P_(C), which evaluates to P_(c)∝½+½ cos [ω(2Δt−δt)]e^(−δt) ² ^(/8σ) ^(coh) ² e^(−Δt) ² ^(/8σ) ^(coh) ² with visibility V=e^(δt) ² ^(/8σ) ^(cor) ² e^(−Δt) ² ^(/8σ) ^(coh) ² . The correlation time and coherence time can be deduced from two visibility measurements V₁ and V₂ using δt₁ and δt₂, respectively. These extrapolated values σ_(coh) ^(E′) and σ_(cor) ^(E′), are given by

$\begin{matrix} {\left( \sigma_{cor}^{E^{\prime}} \right)^{2} = {\frac{1}{8}\frac{{\delta\; t_{1}^{2}} - {\delta\; t_{2}^{2}}}{{\ln\; V_{1}} - {\ln\; V_{2}}}}} & (3) \\ {\left( \sigma_{coh}^{E^{\prime}} \right)^{2} = {\frac{1}{8}\frac{\Delta\;{t^{2}\left( {{\delta\; t_{1}^{2}} + {\delta\; t_{2}^{2}}} \right)}}{{\delta\; t_{1}^{2}\ln\; V_{2}} - {\delta\; t_{2}^{2}\ln\; V_{1}}}}} & (4) \end{matrix}$ The parameters

δt

and

δt₂

are set so that

δt

−

δt₂

≈σ_(cor) and |

δt

≠|

δt₂

| to sample the Franson curve at different points.

FIGS. 2( a) and 2(b) illustrate the choice of

δt₁

=0 and

δt₂

≈σ_(cor). Using (σ_(coh) ^(E))²=1/[(σ_(coh) ^(E′))⁻²−σ_(coh) ⁻²] and (σ_(cor) ^(E))²=1/[(σ_(cor) ^(E′))⁻²−σ_(cor) ⁻²] derived from this measurement, the bound on Eve's information per photon is I_(E)≦log₂(σ_(coh)/σ_(coh) ^(E))+log₂(σ_(cor) ^(E)/σ_(cor)), which is the sum of her information obtained from temporal and frequency measurements. While all change in the wave function can be due to Eve's action, the analysis does not assume that Alice and Bob have perfect detectors or receive unaltered photon pairs. A more pessimistic calculation is required for the mutual information between Alice and Bob.

The probability of Alice and Bob measuring in frequency channels i and j and time bins k and l, is p^(ijkl)=|

ζ_(A) ^(i),ζ_(B) ^(j),σ_(bin,A) ^(k),σ_(bin,B) ^(l)|Ψ

|²=|G^(ijkl)|². These coefficients form a joint probability density function which can be used to compute the Shannon entropy S=Σ_(ijkl)p^(ijkl) log p^(ijkl) and mutual information I (A, B)=S(A)+S(B)−S(A, B).

$\begin{matrix} {\mspace{20mu}{{{I\left( {A,B} \right)} = {{{- 2}{\sum\limits_{i}{{\Gamma^{ik}}^{2}\log{\Gamma^{ik}}^{2}}}} + {\sum\limits_{i,j}{{G^{ijkl}}^{2}\log{G^{ijkl}}^{2}}}}}{\Gamma_{ik} = {\int_{k\;\sigma_{{bi}\; n}}^{{({k + 1})}\sigma_{{bi}\; n}}{\int_{- \infty}^{\infty}{{{FT}_{2}\left\lbrack {\int_{i\;\delta\; v}^{{({i + 1})}\delta\; v}{\int_{- \infty}^{\infty}{{\psi\left( {\omega_{A},\omega_{B}} \right)}{\mathbb{d}\omega_{A}}{\mathbb{d}\omega_{B}}}}} \right\rbrack}{\mathbb{d}t_{A}}{\mathbb{d}t_{B}}}}}}}} & (5) \end{matrix}$

Detector jitter is included by operating a Gaussian spreading function on the biphoton state as {circumflex over (σ)}_(det)=∫e^(−t) ^(x) ² ^(/2σ) ^(det) ² |t

t+t_(x)|dt_(x). FIG. 3 shows the mutual information as a function of the detector jitter normalized to the time bin duration, as discussed hereinbelow. The joint probability density function P of the new state, |Ψ′

={circumflex over (σ)}_(det,A){circumflex over (σ)}_(det,B){Ê}Ŝ_(B)Ŝ_(A)|Ψ_(F)

, incorporating the set of Eve's POVMs, {Ê}, has elements p^(ijkl)=|

ζ_(A) ^(k)ζ_(B) ^(l),σ_(bin,A) ^(i),σ_(bin,B) ^(j)|Ψ′

|².

FIG. 3( a) illustrates mutual information as a function of detector jitter, σ_(det), normalized to σ_(bin). FIG. 3( b) illustrates mutual information as a function of frequency channels, n_(f). The individual frequency channels have bandwidth δv=Δω/n_(f) and spacing 4δv. Increasing the number of frequency channels increases the bits per photon, until the filter photon temporal bandwidth approaches σ_(bin).

The number of time- and frequency-encoded bits given a photon budget and channel bandwidth can be optimized. For example, the photon budget is specific to the photon pair source and is given by the maximum emission rate, R_(v). The allotted channel frequency bandwidth is ΔΩ. The frequency bandwidth of an individual photon pair source is Δω∝1/σ_(cor), so the number of sources to fill the channel bandwidth is N=ΔΩ/Δω. The communication rate, R_(C), in bits per second is therefore

$\begin{matrix} {R_{C} = {{NR}_{v}{\log_{2}\left( \frac{\Delta\;{\omega/\delta}\; v}{R_{v}\sigma_{bin}} \right)}}} & (6) \end{matrix}$ using σ_(bin) and δv defined above.

The often limited photon budget for quantum communication makes high-dimensional encoding desirable. However, achieving the limit on this dimensionality in the time domain using energy-time entangled photon pairs requires detectors with sub-ps timing jitter and resolution. By invoking conjugate frequency correlations, a protocol to approach this fundamental limit has been developed using current detectors and existing telecom networks. The conjugate nature of time and energy encoding means that one can trade frequency for temporal bits (and vice versa) to minimize the effect of channel distortion such as nonlinear frequency conversion and dispersion, in addition to optimizing over transmission rate and channel bandwidth.

Gaussian filtering functions are assumed with center frequencies v_(i) and bandwidth δv, which project onto spatial modes |ζ_(i)

. This is summarized in the following relation, where the total filtering operator Ŝ is taken as a sum over the individual channel filters Ŝ^(i).

$\begin{matrix} {\hat{S} = {{\sum\limits_{i}{\hat{S}}^{i}} = {\sum\limits_{i}{\int_{- \infty}^{\infty}{{\exp\left\lbrack \frac{- \left( {\omega - v_{i}} \right)^{2}}{2\delta\; v^{2}} \right\rbrack}\left. {\omega,\zeta^{i}} \right\rangle\left\langle {\omega,o} \right.{\mathbb{d}\omega}}}}}} & (7) \end{matrix}$

The Franson interference derived hereinabove assumes lossless propagation through the interferometer. However, this assumption is not valid in photonic integrated chips or fiber networks. Loss can be taken into account by placing a beam splitter in the long path of the Franson, which couples the waveguide mode with a vacuum mode (see FIG. 4). Working in the Heisenberg construction, evolving the annihilation operator through the virtual loss beam splitter and the two Franson beam splitters. The matrix for beam splitters 1 and 2, which leave the third mode undisturbed is given by

$\begin{matrix} {{\hat{U}}_{i} = \begin{pmatrix} \sqrt{r_{i}} & \sqrt{1 - r_{i}} & 0 \\ \sqrt{1 - r_{i}} & {- \sqrt{r_{i}}} & 0 \\ 0 & 0 & 1 \end{pmatrix}} & (8) \end{matrix}$ where i ε 1,2. The virtual loss beam splitter is given by

$\begin{matrix} {{\hat{U}}_{L} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & \sqrt{r_{L}} & \sqrt{1 - r_{L}} \\ 0 & \sqrt{1 - r_{L}} & {- \sqrt{r_{L}}} \end{pmatrix}} & (9) \end{matrix}$

The resulting annihilation operators are then â_(A)(t_(A))=C₁â(t)+C₂â(t−Δt) and â_(B)(t_(B))=C₁â(t)+C₂â(t−Δt−δt), disregarding the vacuum term, which will not affect coincidence counting. C₁=√{square root over (r₁)}√{square root over (r₂)} and C₂=√{square root over (1−r₁)}√{square root over (1−r₂)}√{square root over (r_(L))}. For maximum visibility, C₁=C₂, so

$\begin{matrix} {\frac{\sqrt{r_{1}}\sqrt{r_{2\;}}}{\sqrt{1 - r_{1}}\sqrt{1 - r_{2}}} = \sqrt{r_{L}}} & (10) \end{matrix}$

This is plotted. For, r₁=r₂=½, the visibility simplifies to

$\begin{matrix} {V_{PIC} = {\frac{2{\mathbb{e}}^{{- 2}\Delta\;{t/\tau_{\alpha}}}}{1 + {\mathbb{e}}^{{- 4}\Delta\;{t/\tau_{\alpha}}}}{\mathbb{e}}^{{- \delta}\;{t^{2}/2}\sigma_{cor}^{2}}{\mathbb{e}}^{{- \Delta}\;{t^{2}/2}\sigma_{coh}^{2}}}} & (11) \end{matrix}$ where τ_(α) is the photon lifetime in the waveguide due to loss.

Detector jitter refers to the added uncertainty in the photon detection time of some stimulus, purely a result of detector electronics. Superconducting nanowire single photon detectors and InGaAs APDs both exhibit jitter of roughly 30 to 40 ps. Detector jitter is modeled as a Gaussian projection, {circumflex over (σ)}_(det)=∫e^(−t) ^(x) ² ^(/2σ) ^(det) ² |t

t+t_(x)|dt_(z). The jitter profile is not truly Gaussian and can be quite asymmetric, however this model allows for first-order analysis. If this is operated on both Alice and Bob's photons, assuming Eq. (1),

$\begin{matrix} {{{\hat{\sigma}}_{\det,A}{\hat{\sigma}}_{\det,B}\left. \Psi \right\rangle} \propto {\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{{\exp\left\lbrack \frac{- \left( {t_{A} + t_{B}} \right)^{2}}{{4\sigma_{\det}^{2}} + {16\sigma_{coh}^{2}}} \right\rbrack}{\exp\left\lbrack \frac{- \left( {t_{A} + t_{B}} \right)^{2}}{{4\sigma_{\det}^{2}} + {4\sigma_{cor}^{2}}} \right\rbrack}\left. {t_{A},t_{B},o_{A},o_{B}} \right\rangle{\mathbb{d}t_{A}}{\mathbb{d}t_{B}}}}}} & (12) \end{matrix}$

Since σ_(coh)

σ_(det), the most important effect of jitter is to increase the observed correlation time roughly from σ_(cor) to σ_(det). This can have a significant effect on the mutual information between Alice and Bob if σ_(det) is on the order of σ_(bin), as shown in FIG. 3.

The case of a single eavesdropper measuring a single photon of the photon pair is considered herein. Eve's temporal measurement is a POVM that can be written as a Gaussian filtering function: {circumflex over (E)}_(t)=∫_(−∞) ^(∞) e ^(−t) ² ^(/2(σ) ^(coh) ^(E′) ⁾ ² |t

t|dt  (13)

Following, the amplitude function ψ(t_(A),t_(B))αexp[−(t_(A)−t_(B))²/4σ_(cor) ²]exp[−t_(A) ²/4σ_(coh) ²] for σ_(coh)>>σ_(cor). Therefore

$\begin{matrix} {\left. \Psi_{E} \right\rangle = {{{\hat{E}}_{t}\left. \Psi \right\rangle} \propto {\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{{\exp\left\lbrack {- {t_{A}^{2}\left( {\frac{1}{4\sigma_{coh}^{2}} + \frac{1}{4\left( \sigma_{coh}^{E^{\prime}} \right)^{2\;}}} \right)}} \right\rbrack}{\exp\left\lbrack \frac{- \left( {t_{A} - t_{B}} \right)^{2}}{4\sigma_{cor}^{2}} \right\rbrack}\left. {t_{A},t_{B},o_{A},o_{B}} \right\rangle{\mathbb{d}t_{A}}{\mathbb{d}t_{B}}}}}}} & (14) \end{matrix}$ so the coherence time of the biphoton packet is strongly influenced by Eve's filtering bandwidth for σ_(coh) ^(E)>>σ_(coh), which gives the bound on her timing information.

Similarly, a weak frequency POVM is defined, Ê _(ω)=∫_(−∞) ^(∞) e ^(−(σ) ^(cor) ^(E′) ⁾ ² ^((ω−ω) ^(p) ^(/2)) ² |ω

ω|dω  (15)

For 1/σ_(cor)>>1/σ_(coh), and ignoring the spatial modes, |Ψ_(F)

can be written as |Ψ_(F)

≈∫∫exp[−σ_(cor) ²/4(2ω_(A)−ω_(p))²]exp[−σ_(coh) ²(ω_(A)+ω_(B)−ω_(p))²]|ω_(A),ω_(B)

dω _(A) dω _(B).  (16) {circumflex over (E)}_(ω)|Ψ_(P)

≈∫∫exp[−(σ_(cor) ²/4+(σ_(cor) ^(E′))²/4)(2ω_(A)−ω_(p))²]×exp[−σ_(coh) ²(ω_(A)+ω_(B)−ω_(p))²]|ω_(A), ω_(B)

dω _(A) dω _(B).  (17) Thus, Eve projects the biphoton pair onto a narrower frequency distribution.

As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosed subject matter. 

What is claimed is:
 1. A method for distributing a quantum key, the method comprising: sending a first photon to a first receiver; sending a second photon to a second receiver, the first and second photons being a pair of time-energy entangled photons; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins, wherein a combination of a time bin and a frequency bin corresponds to a character.
 2. The method of claim 1 wherein the first photon and the second photon are sent via optical fiber.
 3. The method of claim 1 wherein the first photon and the second photon are sent via a photonic integrated chip.
 4. The method of claim 1 wherein the first photon and the second photon are sent through free space.
 5. A method for distributing a quantum key, comprising: generating a plurality of pairs of time-energy entangled photons; for each pair of time-energy entangled photons, sending one photon to the first receiver and one photon to the second receiver; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins wherein each pair of one of the plurality of time bins and one of the plurality of frequency bins corresponds to one of a plurality of characters.
 6. A method of receiving a quantum key, comprising: receiving a first photon, the first photon being one of a pair of time-energy entangled photons; detecting a time of arrival of the first photon; detecting a frequency of the first photon; and determining a character based on the detected time and frequency of the first photon.
 7. The method of claim 6 wherein the determining a character step further comprises: assigning the detected time to a time bin; assigning the detected frequency to a frequency bin; and determining the character from a coding scheme based on the time bin and the frequency bin.
 8. The method of claim 6, further comprising: detecting whether an eavesdropper has detected the first photon by determining the frequency distribution of the first photon.
 9. A method of receiving a quantum key, comprising: receiving a plurality of photons, each of the plurality of photons being one of a pair of time-energy entangled photons; detecting the time of arrival of each of the plurality of photons; detecting the frequency of each of the plurality of photons; and determining a character based on the detected time and frequency of each of the plurality of photons. 